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For applications to sensor design, the product nx/i of the electron density n and the mobility /i 
is a key parameter to be optimized for enhanced device sensitivity. We model the carrier mobility 
in a two dimensional electron gas (2DEG) layer developed in a ^-doped heterostructure. The sub- 
band energy levels, electron wave functions, and the band-edge profile are obtained by numerically 
solving the Schrodinger and Poisson equations self- consistently. The electron mobility is calculated 
by including contributions of scattering from ionized impurities, the background neutral impurities, 
the deformation potential acoustic phonons, and the polar optical phonons. We calculate the de- 
pendencies of nx/i on temperature, spacer layer thickness, doping density, and the quantum well 
thickness. The model is applied to ^-doped quantum well heterostructures of AllnSb-InSb. At low 
temperature, mobilities as high as 1.3xlO'^m^/V-s are calculated for large spacer layers (400 A) and 
well widths (400 A). The corresponding room temperature mobility is 10 m^/V-s. The dependence 

nxfi shows a maximum for a spacer thickness of 300 A for higher background impurity densi- 
ties while it continues to increase monotonically for lower background impurity densities; this has 
implications for sensor design. 
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I. INTRODUCTION 

In a number of magnetic sensor applications [1-2] the 
basic considerations for improved sensor performance are 
using ultra-thin films, of thicknesses < 100 nm, and very 
high (> 1 m^/V-s) room temperature mobility. However, 
it is also important to optimize the product of the carrier 
density, n, and the mobility, /i, in these devices since they 
have sensitivity that is proportional to the square of the 
mobility [3-4], and because high electron concentrations 
n and electron mobility /igend to reduce the temperature 
variation of the sensor output [1, 4]. 

The mobility in the lateral or in-plane direction, 
in a layered heterostructure, is substantially enhanced 
through selectively doping specific regions [5]. Consider 
a quantum well made up of two layers of compound semi- 
conductor alloys, acting as the barrier regions, on either 
side of a layer with a smaller energy bandgap. On dop- 
ing the barrier regions with impurities, the bound car- 
riers are released into the quantum well which has ac- 
cessible energy levels that are lower than the impurity 
donor energy levels in the barrier region. If the barrier 
region has impurity atoms only on one side, the released 
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free carriers are confined to a narrow region close to the 
barrier- well interface, forming a two-dimensional electron 
gas (2DEG) [6]. The local separation of the charges leads 
to substantial bending of the conduction band-edge, as 
in figure 1. The carriers occupy sub-band energy levels 
in the quantum well and are held in the layer by the pos- 
itively charged ionized donors left behind in the barrier. 
The ionized impurity scattering in lateral transport in 
the quantum well layer, due to the placing of the impuri- 
ties in the barrier regions, can be substantially decreased 
further by set-back doping [7]. 

In this paper, we are concerned with calculating the 
dependence of the mobility on the spacer layer thick- 
ness ds, quantum well thickness d^^ temperature T and 
the doping density in a delta-doped layer in the bar- 
rier. The use of delta-doping with a spacer layer leads 
to a very large mobility for the carriers. We also work 
with the AllnSb-InSb system in which the band gap and 
the electron effective mass in the lateral direction are 
small, which is an advantage for obtaining enhanced mo- 
bility as high as 4.1 m^/V-s at room temperature [8], 20.9 
mVV-s at 77 K [8], and 28.0 mVV-s at 7 K [9], in ul- 
tra thin-film structures. This attractive property based 
on the presence of the InSb layer has been used in ob- 
taining extraordinarily high magnetoresistance (EMR) in 
semiconductor- metal hybrid structures [4]. While a num- 
ber of papers in the literature treat the AlGaAs-GaAs 
[10-14] system theoretically, there are only a few that deal 
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FIG. 1: Figure 1. The conduction band bending, the ground 
state energy level and its wave function in the S -doped layered 
heterostructure are shown. 
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FIG. 2: Figure 2. A schematic diagram of a ^-doped 
Alo . 85 Ino . 15 Sb-InSb heterostructure . 



with transport properties of the AllnSb-InSb structures 
[15]. 

In the following, in Sec. II, we describe the calcula- 
tion of the energy levels in the heterostructure by solving 
Schrodinger's equation, and the corresponding wavefunc- 
tions, together with the solution of the Poisson equation 
for the band edge potential given the redistribution of 
the charges. In Sec. Ill, we evaluate the carrier mobil- 
ity by determining the lifetimes for scattering from the 
impurities and from phonons. Concluding remarks are 
relegated to Sec. IV. 



charge density together with the ionized donor charge 
density at the delta-doped layer are used as the source 
terms for Poisson's equation. The change in the band 
edge profile due to the redistribution of the charges is 
solved for and a fraction of the change (under-relaxation) 
in the potential is added to the Schrodinger's equation to 
initiate the next iteration. At each iteration, the Fermi 
level is determined through the charge neutrality of the 
structure as a whole. Convergence is reached when the 
potential function V changes by less than 0.1 meV over 
the last iteration [12]. 

2.1 Solution of the Schrodinger and the Poisson 
equations 

We consider the AllnSb-InSb structure shown in figure 
2. The donors in the (5-doped layer release electrons into 
the quantum well layer. These electrons are held close to 
the quantum well interface by the ionized impurities left 
behind and form a 2DEG in the triangle-shaped potential 
created by this redistribution of the charges. 

In the effective mass approximation [16-17] we assume 
that the envelope functions of the electrons are propagat- 
ing solutions in the in-plane (x, 7/)-direction and the z- 
dependent envelope function satisfies the one-band equa- 
tion 
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(2.1) 

Here, mjis the effective mass in each layer, Ecb is the 
original band-edge profile while the potential energy due 
to the redistribution of the charges is V{z). We neglect 
nonparabolicity [18] of the conduction band energy dis- 
persion. We solve the above equation by the finite ele- 
ment method (FEM) [19-20]. 

The one dimensional Poisson equation is 
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where Er is the relative dielectric constant. is the ion- 
ized donor concentration, 71(2;) is the electron density in 
the Hartree approximation, and pe is the charge distribu- 
tion [12, 20]. The electron density distribution is related 
to the wave function and subband electron occupation by 



n{z) = J^n^\^^{z) 



(2.3) 



II. SOLVING THE SCHRODINGER AND 
POISSON EQUATIONS 

We first provide an outline of the method. The 
Schrodinger and Poisson equations are solved iteratively 
by first starting with the original band edge profile with 
no band bending. The wavefunctions obtained for the 
electron in the original potential are used to determine 
the charge distribution of the electrons, and this carrier 
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where ^/ is the Fermi energy and Ei is the subband en- 
ergy of the bound states in the quantum well, and B is 
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the step- function. The self-consistent potential V{z) de- 
pends on pe which in turn depends on Ei and tpi obtained 
from the solution of the Schrodinger equation. 

We seek bound state solutions to the Schrodinger equa- 
tion with V^(0) = at z= at the top of the het- 
erostructure and ^Ij{L) = Oat the bottom. At layer 
interfaces, the wavefunction continuity and the mass- 
derivative of the wavefunction, tjj' /m^^ are assumed to 
be continuous [21]. The boundary conditions for the Pois- 
son equation are: V{z =0) = Vc, a contact poten- 
tial, and V{z = L) = 0, at the bottom of the het- 
erostructure. At layer interfaces, we require that V and 
SrdV/dz be continuous. In order to solve Eqs. (2.1- 
2.4) self-consistently, we start with the initial choice for 
the potential function, V= 0, and solve the Schrodinger 
equation for the bound state energies £n and wavefunc- 
tions ^/^^^^.From this we construct the first approxima- 
tion to the charge densitype- Now Poisson's equation is 
solved to determine the change in the potential energy, 
V(z), with the source terms given by the charge den- 
sity — e [n {z) — Nd {z)]. A fraction of this potential en- 
ergy function is added to the potential energy used in the 
Schrodinger equation. Letting k be the iteration index, 
we set F^+i {z)=a' {z) + (1 - a) • V^'^ (z), where a 
is ~0.01. This is used, in turn, in the Schrodinger equa- 
tion for the next iteration of the wavefunction if)^^ using 



(2.5) 

These wavefunctions are used in defining n^^^(yz) in the 
source term in Poisson's equation 
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The process is repeated until convergence is reached 
through the criterion that the change in the potential 
function reaches the tolerance A such that 



|V^+i {z)-V^ {z)\ 



< A. 



(2.7) 



The tolerance was chosen to be A = 1 x 10~^. The re- 
sulting conduction band edge profile is shown in figure 1, 
together with the ground state wavefunction. The con- 
duction band edge Ecb in each layer is determined by a 
parameter, 7, representing the conduction band offset in 
the layer with respect to the band-gap, Eg^ in the ma- 
terial and a reference layer. We include the temperature 
dependence of the band gap as given by the empirical 
Varshni equation [22] 



E, (T) = Eg{T = 0) 



T + 13 



(2.8) 



where a and (3 are the usual Varshni parameters [23]. 



For all the ternary alloys, the composition-dependent 
bandgap is given by a simple linear or quadratic equation 
as a function of concentration x by [24-25] 



Eg {Ai_,B,) = {l-x) Eg {A) + xEg {B) -x{l 



■x) C. 
(2.9) 

Where C is the bowing parameter for the bandgap, and 
accounts for the deviation from a linear interpolation 
between the bandgaps of the two binaries A and 5 [23]. 

2. 2 The determination of the Fermi level 

As the electrons fall into the quantum well the band- 
edge potential in the well region moves upwards to make 
it less energetically unfavorable for the next electron to 
enter the well region. If band nonparabolicity is ignored, 
the 2DEG density is given by [11] 
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If the barrier doping density is large, the location of the 
Fermi level is defined by the energy level of the unionized 
donor and charge neutrality requires that 



(2.11) 



The Fermi level can be determined by using a root-finding 
procedure. Equation (2.11) is solved for the Fermi level 
in every iteration, in order to speed up the convergence. 



III. CALCULATION OF THE MOBILITY 

The carrier mobility, /i = er/m*, is governed by 
scattering mechanisms that control the scattering time, 
r, associated with each scattering process in a 2DEG 
[26] . The scattering rates calculated with each scattering 
mechanism are combined to determine the resultant 
mobility. The dominant scattering mechanisms for bulk 
III-V compounds are now well established [27-29]. In 
our calculations, we have included these mechanisms in 
the context of 2D scattering. 

3. 1 The scattering from ionized impurities 

At low temperature, scattering is dominated by the 
ionized donors. In a (5-doped structure with a spacer 
layer, the ionized donors and the 2D electrons are spa- 
tially separated, thereby minimizing the scattering. We 
assume that the electron gas is highly degenerate and 
that the scattering occurs only with electrons near the 
Fermi level. The transport scattering rate for a purely 
2D electron gas, neglecting the electronic wavefunction 
normal to the plane of the 2D gas, is given by [16] 
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where kp = ^Jn\^ • 27r is the Fermi wavevector, v?f^ is Fourier transform of the scattering potential. 

the sheet density of the 2DEG, A/',^il, is the impurity dop- • • • ^ r 

^ ' ^"^^ ^ Usmg ihomas-l^ermi screenmg, the scattermg rate tor 

ing areal density and V (^) = (r) e~'^ ' ^ (P r is the remote impurities becomes [16] 
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and b = (^^i^^^^) • The scattering potential takes 
account of screening with the Thomas-Fermi dielectric 

function [30]. Here qrr = ^txK^ee screening 
wavevector in the Thomas-Fermi approximation. The 
contribution to the mobility from ionized impurity 
scattering is essentially independent of temperature [31]. 



S.2 The ionized background impurities 

The background density of ionized donor impurities is 
usually very small as compared with ionized impurities 
with thin spacer layer thickness. But when the thickness 
of the spacer layer is very large, the effect of remote ion- 
ized impurities is reduced and the background impurity 
scattering dominates. The rate for scattering from the 
low density background impurities with a concentration 
ofA^5'™ is given by [16] 



S. S The scattering from phonons 

Phonons dominate scattering at high temperature, 
typically above 60 K [32]. The most important phonon 
scattering processes in general are: {i) deformation po- 
tential acoustic phonon scattering, (ii) polar optical 
phonon scattering, and (Hi) piezoelectric scattering [33- 
34]. However, Basu and Nag [35] have shown that for 
InSb, the piezoelectric scattering does not play an impor- 
tant role at intermediate temperatures. We will therefore 
not consider it here. We consider the electrons to be 
quasi-two dimensional while the phonons remain quasi- 
three-dimensional, in which approximation the perturb- 
ing potential has a spherical symmetry [16]. For the sake 
of simplicity, we assume that in the ideal case we have at 
hand bulk phonons. 

The deformation potential acoustic phonon scattering 
rate is given by [32, 36] 
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where S is the deformation potential, and p, Vs are the 



mass density and the velocity of sound. 

The scattering rate by polar optic phonons is given by 
[32-33] 



koe'^nNo 



polar 



(3.5) 



with N( 







(e^-l) , ^0 = V^^, and 
— = — — — , with being the high frequency 
dielectric constant and £r being the static dielectric 
constant. Here hujo is the optical-phonon energy. Several 
values (7.2 eV and 30 eV [37-38]) have been quoted 
in the literature for the deformation potential in InSb. 
Since the mobility is inversely proportional to S^, these 
different S values can lead to a large difference in the 
acoustic phonon scattering contribution to the mobility. 
We use a value of 7.2 eV for S in InSb, which was first 
obtained by Ehrenreich [39] and Dutta [38] who showed 
that it agrees well with the experimental results. 

3.4 Other scattering mechanisms 
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The alloy- disorder scattering by free carriers was 
neglected in this calculation. Interface (or surface) 
roughness scattering is another scattering mechanism 
which has been found to be important in thin quantum 
wells [40], playing a significant role only at high electron 
concentrations [41], and thin quantum wells ( L<60 A) 
[42] . When the first excited subbands become filled with 
electrons, there will be interface scattering, which would 
reduce the mobility [43]. 

3.5 The combined mobility 

The total mobility is given by computing the mobil- 
ity for each scattering process and adding the reciprocal 
mobility for each process according to Mathieson's rule. 



1 



1 



l^total l^imp 



1 



1 



f^acoustic f-^polar 



(3.6) 



At high temperature, the relation may not be valid due 
to the limited applicability of degenerate statistics [31]. 
However, for temperatures higher than 60 K the scatter- 
ing in the 2DEG is dominated by the phonon scattering 
processes. The approximation has a negligible effect on 
the combined electron mobility. Therefore we still can 
use Mathieson's rule to roughly estimate the combined 
mobility. 

The mobility contributions caused by ionized impurity 
and background impurity scatterings at low temperatures 
are temperature independent, above the temperature at 
which the total mobility shows the T~^ dependence. 



IV. 



RESULTS AND DISCUSSION 



We studied the Alo.SSIno.isSb-InSb system. Some 
relevant parameters used in the calculation are given in 
Table 1. The other material parameters employed here 
have been obtained from Ref. 23. Table 1. Parameters 
used in the calculation are shown. 



Parameters of InSb 


Value 


Effective mass ratio (m*/mo) 


0.0135 


Static dielectric constant {so ) 


16.82 


High-frequency dielectric constant {soo ) 


15.7 


Lattice density (g/cm^) 


5.79 


Velocity of longitudinal elastic waves (cm/s) 


3.7x10^ 


Acoustic deformation potential (eV) 


7.2 


Optical phonon energy (meV) 


25.0 



^) Ref. 23 ^) Ref. 15 Ref. 37 

4.1 Temperature dependence of the electron mobility 

The calculated electron mobility, in the range of 
1-280 K, is given for the AllnSb-InSb system in figure 
3. The components of the mobility contributing to 
the total mobility are also presented. In figure 3(a), 
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FIG. 3: Figure 3. The temperature dependence of the electron 
mobility in the Alo.ssIno.isSb-InSb system is shown, (a) The 
curves show the calculated mobility with scattering due to 
the remote ionized impurities with a sheet density of 1 x 10^^ 
cm^. The spacer layer thickness and well thickness are 400 A. 
The background impurity density is assumed to be 5x10^^ 
cm~^. (b) The temperature dependence of the electron sub- 
band density and the product of the carrier density and the 
mobility, nx/i, are shown. 



with T decreasing from 280 K, the carrier mobility 
increases monotonically and saturates around 30 K. 
This is in contrast to the bulk semiconductor, where the 
mobilities tend to peak at intermediated temperature. 
This is due to the fact that at low temperature, the 
electron mobility is limited by ionized impurity scatter- 
ing and phonon-limited scattering is negligible; on the 
other hand, at high temperature polar-optical-phonon 
scattering is the dominant scattering mechanism. At 
intermediate temperature, the deformation-potential 
acoustic phonon scattering also plays an important role. 
At high temperature, the mobility shows the expected 
T~^ dependence. The subband electron occupation 
is shown in figure 3(b). In our case, only the ground 
state subband is filled. The electron density in the 
quantum well remains constant at low temperature. But 
it increases rapidly about 100 K. We know that for the 
sensor application, high electron density determines the 
conductivity and this is as important as the need for 
high electron mobility. We have therefore plotted the 
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FIG. 4: Figure 4. (a)The dependence of the electron mobiUty 
on the thickness of the spacer layer in the Alo.ssIno.isSb-InSb 
heterostructure is shown. The calculation is performed at 
T=0 K for a well width of 200 A. The background impurity 
density is assumed to be SxlO^^cm""^. (b) The dependence 
of electron subband density and the product of the carrier 
density and the mobility, nx/i, on the spacer layer thickness 
are shown. 
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FIG. 5: Figure 5. (a) The dependence of the electron mo- 
bility on the spacer layer thickness in the A^o.ssIno.isSb-InSb 
heterostructure is shown. The background impurity density 
is assumed to be 5xl0^^cm~^. (b) The product of the carrier 
density and the mobility, nx/x, is shown as a function of the 
spacer layer thickness. 



product of electron density n and mobility /i in figure 
3(b). The product of nx/i decreases rapidly at high 
temperature, and remains constant at low temperature. 

4-2 Spacer thickness dependence 

The purpose of a (5-doped heterostructure is to separate 
the 2DEG from the parent ionized donors, thereby lim- 
iting ionized impurity scattering from the remote doping 
centers. Figure 4(a) shows that an increase in the spacer 
thickness leads to an increase in the combined mobility. 
The ionized impurity mobility increases with increasing 
spacer thickness, reaching mobilities of 1x10^ m^/V-s at 
a spacer width of 400 A. The mobility caused by the in- 
troduction of background ionized impurities shows the 
reverse trend. The background impurity mobility de- 
creases with increasing spacer width and becomes domi- 
nant. However, an increase of spacer layer thickness leads 
to a decrease in the electron density n in the quantum 
well as shown in figure 4(b). The product nx/i and the 
carrier density n are also shown in figure 4(b), with the 
former increasing with spacer layer thickness. 

For a background density of 5x10^^ cm~^, the spacer 



width dependence of combined mobility is insensitive to 
the background impurity as can be seen in figure 4(a). 
However, higher background impurity concentrations 
will make the background impurity mobility comparable 
with the ionized impurity mobility which will result in 
the combined mobility displaying a peak at some spacer 
thickness, as shown in figure 5. Figure 5(b) also shows 
that nx/i as a peak value at a spacer width around 300 
A at high background impurity density. 

4.3 Dependence on the quantum well width 

For the well width dependence study we considered the 
system with a constant spacer layer thickness of 200 A. 
Figures 6(a) and (b) show that at K for a thin quan- 
tum well ( < 200 A), electron density nand nx/i increase 
rapidly up to 200 A, while at larger well width, 

they increase slowly and saturate at a well width 300 A. 
With the increase in the well width, the background im- 
purity mobility increases faster than the ionized impurity 
mobility and is a less important factor in the combined 
mobility. The electron density n and the product nx/i 
show the same trend with increasing well width. 

Equations (3.4)- (3.6) show that the acoustic-phonon 
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FIG. 6: Figure 6. (a) The dependence of the electron mobihty 
on the well width in the Alo.ssIno.isSb-InSb heterostructure is 
shown. The calculation is performed at T=0 K and a spacer 
layer thickness of 200 A. A background impurity density of 
5x10^^ cm~^ is assumed, (b) The dependence of electron 
subband density, and of the product of the mobility and elec- 
trons, nx/i, on the well layer thickness is shown. 
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FIG. 7: Figure 7. The dependence of the electron mobility on 
the well layer thickness in the Alo.ssIno.isSb-InSb heterostruc- 
ture is shown on a log- log plot, so that the dependence of 
polar phonon scattering on well width is shown to be nearly 
constant. The calculation is performed with T=40 K and a 
spacer thickness of 300 A. 
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FIG. 8: Figure 8. (a) The dependence of the electron mobility 
on the doping density in the Alo.ssIno.isSb-InSb heterostruc- 
ture is shown. The calculation is performed with T=0 K and 
spacer and a well thickness of 200 A. The background impu- 
rity density is assumed to be 5xl0^^cm~^. (b) The variation 
of the Fermi level, the ground subband energy and electron 
densities with the ionized doping density are shown, (c) The 
product of mobility and electron density, nx/i, versus the dop- 
ing density is shown. 



mobility should increase linearly with the well width 
at intermediate temperature, and that the polar optic 
phonon scattering remains constant. Figure 7 shows the 
calculated mobility at 40 K. Background impurity and 
ionized impurity mobility show a similar trend compared 
to the K mobility (Fig. 6a). 



^.^ Doling density dependence 
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Figure 8 shows that at K as the remote impurity 
density increases, the background impurity mobihty in- 
creases faster than the ionized impurity mobihty. This 
fohows logicahy, because with more electrons in the quan- 
tum weh, the effect of background impurities is not the 
hmiting factor in the combined mobihty, therefore it 
shows a trend similar to that of the electron subband 
density in figure 8(b). At low doping density, all the 
electrons can enter the well without pushing the Fermi 
energy up to the conduction band edge. With increased 
doping density, more electrons can enter the well, and 
the Fermi energy reaches the conduction band-edge and 
is unable to move up further. The ground subband en- 
ergy will then be lower to let more electrons into the well 
as shown in figure 8(b). 

V. CONCLUSION 

In conclusion, we have presented theoretical calcula- 
tions for the electron subband energy, subband electron 
density, and carrier mobility in S -doped single quantum 
well heterostructures. Here, the Hartree approximation 
for the confinement potential was used. The important 
scattering mechanisms, such as, ionized impurity scatter- 
ing, background impurity scattering, deformation poten- 
tial acoustic phonon scattering, and polar phonon scat- 



tering were considered. We have presented the behavior 
of the product of electron mobility and density in the 
temperature range of - 280 K as a function of the spacer 
layer thickness, well width, and doping density. 

At low temperature, the dominant scattering mecha- 
nism is ionized impurity scattering. It is also found to 
be independent of temperature in this regime. In the 
AllnSb-InSb system, the mobility due to ionized impurity 
scattering increases with decreasing temperature, reach- 
ing a limiting value of about 1x10^ cm^/V-s for a thicker 
well and larger spacer ( >300 A) at low temperature ( 
K). At a temperature around 40 K, optical phonons begin 
to limit the mobility which varies as T~^at high temper- 
ature. The dependence of nx/u shows a maximum for a 
spacer thickness of 300 A for higher background impurity 
densities while it continues to increase monotonically for 
lower background impurity densities. 
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